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Aerothermodynamic design environments are presented for the Mars Science Labora- 
tory entry capsule heatshield. The design conditions are based on Navier-Stokes flowfield 
simulations on shallow (maximum total heat load) and steep (maximum heat flux, shear 
stress, and pressure) entry trajectories from a 2009 launch. Boundary layer transition is 
expected prior to peak heat flux, a first for Mars entry, and the heatshield environments 
were defined for a fully-turbulent heat pulse. The effects of distributed surface roughness 
on turbulent heat flux and shear stress peaks are included using empirical correlations. 
Additional biases and uncertainties are based on computational model comparisons with 
experimental data and sensitivity studies. The peak design conditions are 197 W/cm 1 2 for 
heat flux, 471 Pa for shear stress, 0.371 Earth atm for pressure, and 5477 J/cm 2 for total 
heat load. Time- varying conditions at fixed heatshield locations were generated for thermal 
protection system analysis and flight instrumentation development. Finally, the aerother- 
modynamic effects of delaying launch until 2011 are previewed. 


Nomenclature 


A reference area, \ttD 2 (m 2 ) 

Cd drag coefficient, D/q^A 

D aeroshell diameter (m) 

Dim multi-component diffusion coefficient (m 2 / s) 

Ci species mass fraction 

H total enthalpy (J /kg) 

hi species enthalpy (J /kg) 

k roughness height (mm) 

k + roughness Reynolds number, p w U T k s / fi w 

k c thermal conductivity (W/m — K) 

k s equivalent sand-grain roughness height (mm) 

L/D lift-to-drag ratio 

m entry system mass (kg) 

p pressure (Earth atm , 1 Earth atm = 101,325 Pa) 

Q convective total heat load, f qdt (J/cm 2 ) 

q convective heat flux (W/cm 2 ) 

goo freestream dynamic pressure, ^Poo^io {P°) 

Re oo freestream Reynolds number, Poofoo^/Poo 

* Aerospace Engineer, Atmospheric Flight and Entry Systems Branch, Mail Stop 489, Senior Member, 
Karl.T.Edquist@nasa.gov. 

t Aerospace Engineer, Atmospheric Flight and Entry Systems Branch, Mail Stop 489, Member. 

£ Senior Research Scientist, Reacting Flow Environments Branch, Mail Stop 230-2, Associate Fellow. 

§ Senior Research Scientist, Reacting Flow Environments Branch, Mail Stop 230-2, Member. 


1 of 21 


American Institute of Aeronautics and Astronautics 



Reg 

Sc t 

St 

t 

U T 

U 

V 

X,Y,Z 

Xi 

Xfran 

a 

Pm 

7 


T] 

0 

V 

p 


Acronyms 


AEDC 

AWG 

B-L 

CFD 

DPLR 

LAURA 

MEADS 

MEDLI 

MISP 

MSL 

PICA 

RSS 

SLA 

SST 

TPS 


Subscripts 

D 

e 

trim 

w 

w, 0 
w, k 
oo 


boundary layer momentum thickness Reynolds number, p e u e 6 / fi e 
turbulent Schmidt number, v/D^ 

Stanton number, q w / PooKo(#oo - H w ) 

time from atmospheric interface (s) 

shear velocity, yjr w / p w (m/s) 

tangential velocity (m/s) 

velocity relative to atmosphere (km/s) 

coordinates measured from aeroshell nose (m) 

independent parameter 

Baldwin-Lomax turbulence model transition location (m) 

angle- of- attack (deg) 

ballistic coefficient, m/CoA (kg/m 2 ) 

inertial flight path angle (deg) 

emissivity 

surface normal coordinate (m) 
boundary layer momentum thickness (m) 
kinematic viscosity (m 2 /s) 
density (kg/m 3 ) 
shear stress (Pa) 


Arnold Engineering Development Center 
Aeroheating Working Group 
Baldwin-Lomax 

Computational Fluid Dynamics 
Data Parallel Line Relaxation 

Langley Aerothermodynamic Upwind Relaxation Algorithm 

MEDLI Entry Air Data System 

MSL Entry, Descent, and Landing Instrumentation 

MEDLI Integrated Sensor Plugs 

Mars Science Laboratory 

Phenolic Impregnated Carbon Ablator 

Root Sum Square 

Super Lightweight Ablator 

Shear Stress Transport 

Thermal Protection System 


aeroshell diameter 

boundary layer edge (H e = 0.995iLoo) 
aerodynamic trim condition 
surface condition 
smooth surface condition 
rough surface condition 
freestream condition 


I. Introduction 

NASA’s Mars Science Laboratory (MSL) entry system 1 will deliver a large rover (> 900 kg) to the 
surface of Mars within a distance of 10 km from the landing target. The entry system will land the rover 
using an aeroshell and parachute architecture that has successfully delivered payloads on Mars starting with 
the Viking 2 landers and continuing with Pathfinder, 3 the Mars Exploration Rovers, 4 and Phoenix. 5 Prior to 
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supersonic parachute deployment, the MSL entry capsule will fly with a hypersonic lift-to-drag ratio (. L/D ) 
of 0.24 at a trim angle-of-attack (a tr im) of 16 deg using active guidance and thrusters to control lift vector 
direction. Ballast mass will be used to give the capsule the necessary radial center-of-gravity offset for the 
desired a tr im and L/D. MSL launch, until recently scheduled for the Fall of 2009, has been delayed until 
late 2011 with landing in Summer of 2012. 

The interaction between the MSL aeroshell and Martian atmosphere will dissipate more than 99% of 
the entry system’s initial kinetic energy, mostly in the form of heat. MSL’s heatshield thermal protection 
system 6-8 (TPS) will keep the aeroshell interior from these extreme environments. TPS material response 
is primarily a function of convective heat flux via conduction and diffusion (radiation is negligible for MSL): 


Qw — 


7 W' , 

k c ^j + P / j Ldjmdj 



( 1 ) 


The objective of this paper is to summarize the MSL aerot her mo dynamic design process and conditions 
(heat flux, shear stress, pressure, and total heat load) in support of heatshield TPS analysis 6 and testing, 7 as 
well as development of engineering instrumentation. 9 A companion paper 10 covers design environments for 
the entry capsule backshell and parachute cover. MSL’s large aeroshell and high ballistic coefficient (/3 m ) will 
result in boundary layer turbulent transition prior to peak heating (a first for Mars entry), leading to more 
severe aerothermodynamic conditions than any previous Mars entry heatshield. Consequently, significant 
effort has been spent to quantify the aerothermodynamic environments with ground testing and analytical 
predictions. Computational 11-15 and experimental 16-21 aerothermodynamics analyses have been performed 
for various MSL aeroshell configurations and design trajectories. This paper covers the aerothermodynamic 
analysis used to develop the flight heatshield TPS for the original 2009 launch opportunity. 

A. Aeroshell and Heatshield Thermal Protection System 

The MSL aeroshell is a 4.5-meter diameter spherically-blunted 70 -degree half-angle cone forebody with a 
triconic afterbody (Figure la). All previous major NASA Mars capsules used the same forebody shape and 
Lockheed Martin’s Super Lightweight Ablator 22 (SLA-561V) TPS material on the heatshield. The unprece- 
dented magnitude of the MSL environments led to SLA-561V test anomalies 7 and necessitated switching 
to the Phenolic Impregnated Carbon Ablator 23 (PICA). PICA is developed by Fiber Materials Incorpo- 
rated (FMI) and was first used on the Stardust Earth return capsule, 24 but has never been used on a Mars 
heatshield nor in a tiled configuration (Figure lb). PICA was also under consideration for the Orion Crew 
Exploration Vehicle lunar return capsule. A uniform PICA thickness of 1.25 in was established for MSL 
after the material switch and prior to completion of the aerothermodynamic design environments and TPS 
analysis. Therefore, the thickness was determined by the entry system mass constraints rather than through 
thermal analysis in response to a particular heat pulse. 6 Instead of deriving a minimum required thickness, 
the TPS analysis focused on ensuring that 1.25 in is sufficient to withstand the heat pulse with sufficient 
margin and acceptable risk. The PICA sizing models predict significant recession (more than half of the 
initial thickness) in response to the severe environments, another first for a Mars heatshield. 

Numerous tests 8 were conducted to develop PICA for MSL’s purposes and TPS analyses have since con- 
firmed that 1.25 in is sufficient to meet the thermal requirements with acceptable risk. 6 PICA application has 
been completed for the flight heatshield in anticipation of the original 2009 launch date. PICA was installed 
on the heatshield in tile form (Figure lb) due to manufacturing size constraints. Requirements were placed 
on the gap/flow angles and maximum gap length to minimize adverse effects due to the boundary layer. 
A total of 113 tiles (27 unique shapes) was used to cover the heatshield. Adhesive materials were used to 
attach the tiles to the aeroshell structure and to fill the gaps between tiles. The bondline interface must stay 
below temperature limits (225 to 250° C) to maintain bond integrity and aeroshell structure performance. 6 
The effects of protruding gap fillers were analyzed computationally, 25 but were not explicitly included in 
the design environments because of low probability at high heat fluxes and for gaps perpendicular to the 
oncoming flow 7 where bump factors could pose a risk to the TPS. 
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(a) Geometry (b) PICA Heatshield TPS Tiles 


Figure 1. MSL Aeroshell 


B. Aerothermodynamic Design Trajectories 

Two entry trajectories (08-TPS-01a and 08-TPS- 
02) were used to define MSL’s heatshield aerother- 
modynamic conditions (Figure 2 and Table 1). 
The trajectories were designed with a high en- 
try velocity to give upper bound aerothermal 
loads. The trajectories differ only in en- 

try flight path angle (7). The shallow 08- 
TPS-Ola trajectory (entry 7 = -14 deg ) has 
a longer time of flight and the highest inte- 
grated total heat load, which typically requires 
the largest TPS thickness. The steep 08- 
TPS-02 trajectory (entry 7 = -15.5 deg ) gives 
the peak instantaneous conditions (heat flux, 
shear stress, and pressure) for ground test de- 
sign. The capsule will generate enough lift 
to achieve nearly horizontal flight leading up 
to parachute deployment (between 1 and 2 
km/ s). 



Velocity (km/s) 


Figure 2. Aerothermodynamic Design Trajectories for 
the 2009 Launch Opportunity 


Table 1. Aerothermodynamic Design Trajectories for the 2009 Launch Opportunity 



08-TPS-01a / 08-TPS-02 

Launch Date 

09/20/09 

Arrival Date 

07/20/10 

Entry Mass (kg) 

3380 

Entry Voo (km/s) 

5.63 

Entry 7 (deg) 

-14 / -15.5 

Time at Parachute Deployment (s) 

280 / 249 


4 of 21 


American Institute of Aeronautics and Astronautics 



C. MSL Compared to Previous Missions 


Figure 3 and Table 2 compare MSL to past success- 
ful Mars missions, all of which used SLA-561V in 
less severe heatshield conditions than MSL will ex- 
perience. The unprecedented MSL entry system size 
and mass will lead to a much more aggressive atmo- 
spheric entry than has been previously attempted. 
The MSL rover mass alone (850 kg for the design 
trajectories, > 900 kg for the flight hardware) is sim- 
ilar to the total entry system masses for the Mars 
Exploration Rovers (MER). MSL’s landing mass in- 
cludes the Sky Crane 1 descent stage since it will 
reach the ground after releasing the rover. MSL’s 
high ballistic coefficient leads to higher aerother- 
mal and structural loads since velocity is maintained 
deeper into the atmosphere. Furthermore, MSL’s 
large heatshield surface area enhances the likelihood 
of turbulent transition with longer boundary layer 
running lengths. Finally, MSL will fly at a higher 
angle-of-attack than has ever been attempted to 
achieve the desired L/D ; this further exacerbates 
the aerothermal loads by allowing longer running 
lengths for turbulent boundary layer development. 
MSL’s turbulent heating levels led to SLA-561V test 
failures 7 and the thickest Mars heatshield ever built 



Velocity (km/s) 

Figure 3. Mars Entry Trajectory Comparison 


Table 2. Mars Aeroshell and Entry Trajectory Comparison 

| Viking I fc II | Pathfinder | MER A fc B | Phoenix | MSL (Design) 



Diameter (m) 

3.51 

2.65 

2.65 

2.65 

4.5 

Entry Mass (kg) 

981 

585 

840 

573 

3380 

Entry Voo (km/s) 

4.5 

7.6 

5.5 

5.5 

5.6 

Entry 7 (deg) 

-17 

-13.8 

-11.5 

-13 

-14 to -15.5 

Landing Mass (kg) 

603 

360 

539 

364 

850 + 775 

Hypersonic /l m (kg/m 2 ) 

64 

62 

95 

62 

146 

Hypersonic a trirn (deg) 

11 

0 

0 

0 

16 

Hypersonic (L/D) trirn 

0.18 

0 

0 

0 

0.24 

Turbulent at Peak q w l 

No 

No 

No 

No 

Yes 

Peak q w ( W/cm 2 ) 

22 

118 

48 

55 

197 

Peak p w (atm) 

0.10 

0.20 

0.10 

0.08 

0.37 

Heatshield TPS 

SLA-561V 

SLA-561V 

SLA-561V 

SLA-561V 

PICA 


(0.5 in) 

(0.75 in) 

(0.62 in) 

(0.55 in) 

(1.25 in) 
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D. Computational Models 

There exists no ground-based facility that can reproduce the high-temperature non-equilibrium CO 2 — N 2 
thermodynamic conditions that exist at Mars hypersonic flight conditions. Consequently, extensive use 
of Navier-Stokes computational fluid dynamics (CFD) was required to define the MSL aerot her mo dynam- 
ics. Many challenges 26 exist when calculating the aerothermodynamic environments for a hypersonic Mars 
entry capsule. The high temperatures generated behind the bow shock leads to chemical and thermal non- 
equilibrium conditions that require models whose validity at Mars flight conditions often cannot be verified. 
For MSL, verification and validation of the CFD models was accomplished through code-to-code comparisons 
and anchoring to perfect-gas and high-enthalpy test data. 16-21 

The MSL aerothermodynamic environments are based on the Langley Aerothermodynamic Upwind Re- 
laxation Algorithm 27 (LAURA) and Data Parallel Line Relaxation 28 (DPLR) Navier-Stokes codes. Both 
codes have been used for previous Mars and Earth entry capsule flowfield predictions, including Pathfinder, 29 
Phoenix, 30 and Genesis Sample Return Capsule. 31,32 LAURA is also being used 33 to predict MSL’s aerody- 
namic characteristics in the hypersonic and supersonic continuum flow regimes, including RCS interference 
effects. 34 For MSL, turbulent solutions using the Baldwin-Lomax algebraic model 35 were obtained for the 
heatshield based on the expectation of turbulent transition in flight. It will be shown that the Baldwin-Lomax 
model agrees well with fully-turbulent test data. DPLR solutions were used for verification and validation 
of the LAURA results and sensitivity analyses using the Shear Stress Transport 36 (SST) turbulence model. 
The model assumptions for flight calculations were as follows: 

• Turbulent boundary layer: Baldwin-Lomax algebraic model, transition forced at stagnation point 

• Chemical non-equilibrium (CO 2 , CO, N 2 , O 2 , NO, C, N , O, 0.97 and 0.03 freestream mass fractions for 
CO 2 and N 2 , respectively) 

• Thermal non-equilibrium (translational and vibrational temperatures) 

• Radiative-equilibrium wall: q w = eaT 14 , e = 0.89, cr = 5.67 x 10 -8 W/(m 2 — K 4 ) 

• Super-catalytic wall: recombination of CO 2 and N 2 to freestream mass fractions 

• Non-blowing, smooth wall (distributed roughness effects added post-CFD) 

II. Results 

The following sections summarize the results of applying the process in Figure 4 to define MSL aerother- 
modynamic design environments. Given freestream conditions along the design trajectories, CFD solutions 
were obtained using the model assumptions stated previously. Margins were then added to account for CFD 
model corrections and uncertainties. The end results are the design environments for multiple heatshield 
locations for the 08-TPS-01a (maximum total heat load) and 08-TPS-02 (maximum heat flux, shear stress, 
and pressure) entry trajectories. Any changes that affect the entry trajectory, such as entry mass or velocity, 
required a repeat of the design cycle. Due to the timing of the switch to PICA, the TPS sizing analysis 
turned into a exercise to verify that the as-built thickness (1.25 in) was sufficient rather than to determine 
a minimum required thickness for a particular heat pulse. 6 

MSL aeroshell and heatshield TPS development was a joint effort between NASA Langley and Ames 
Research Centers (aerothermodynamic analysis, TPS modeling and testing), NASA Jet Propulsion Lab- 
oratory (JPL, project management, verification and validation), Lockheed Martin Astronautics (aeroshell 
development), and FMI (PICA manufacturer). The aerothermodynamic analyses that supported the PICA 
design was overseen by the MSL Aeroheating Working Group (AWG) consisting of members from NASA, 
JPL, and Lockheed Martin. The AWG served as the venue for discussing the state of the aerothermal design 
environments, interpreting test results, and proposing additional testing and analysis. The aerothermal de- 
sign process and results were reviewed by project management and subject matter experts, both in detailed 
and system-level project reviews. The data to follow reflect the final analysis vetted through the AWG and 
used to size and test the flight heatshield TPS for the original 2009 launch opportunity. A preview is also 
given of environments for the new launch date in 2011. 
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Figure 4. Aerothermodynamic Design Process 


A. Boundary Layer Transition 

MSL’s large aeroshell, high ballistic coefficient, and non-zero trim angle-of- attack is expected to lead to fore- 
body boundary layer transition prior to peak heating. This assessment was made using a simplified method 
based on exceeding a critical laminar boundary layer momentum thickness Reynolds number : 

p e u e 6 
Re e = 

Me 

where the boundary layer momentum thickness is defined as 



Boundary layer edge conditions are defined where the total enthalpy is 99.5% of the freestream value. 

A conservative transition criterion of Reg > 200 
was established for the heatshield design environ- 
ments. Experimental MSL data in perfect gas and 
high-enthalpy facilities support the use of this tran- 
sition criterion. 20 It is acknowledged that using 
ground-based results to support transition predic- 
tion at flight conditions may be inappropriate. How- 
ever, conservatism was warranted given the uncer- 
tainty in predicting transition onset and its signif- 
icant effect on both heat flux and shear stress for 
a flagship NASA mission. Figure 5 shows that the 
LAURA Reg on the heatshield leeside ( X < 0) ex- 
ceeds 200 prior to the time of peak heat flux on 
the shallow design trajectory. The highest values of 
Reg are found in the leeside shoulder region where 
the boundary layer running length is longest. This 
will be shown to also correspond the area of high- 
est turbulent heat flux and shear stress. In order to 
conservatively account for turbulence effects on total 
heat load, both design trajectory heat pulses were 
assumed to be entirely turbulent. CFD solutions on 
each trajectory were run using the Baldwin-Lomax 
algebraic model with transition forced at the stag- 
nation point. 



X (m) 

Figure 5. LAURA Laminar Momentum-Thickness 
Reynolds Number Distributions on the 08-TPS-01a 
Trajectory 


( 2 ) 
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B. CFD Compared to Experimental Data 

Turbulent transition prior to peak heating is believed to be a first for a Mars entry. Consequently, there exist 
no flight data to validate CFD turbulence models at Mars conditions. For MSL, turbulence model validation 
was achieved through comparison to test data. Several MSL experimental tests 16-21 have been conducted 
to understand boundary layer transition and the validity of the Baldwin-Lomax model. One of those test 
programs 16 was conducted in the Arnold Engineering Development Center (AEDC) Tunnel 9 using perfect 
gas nitrogen. The facility’s wide Reynolds number range allowed data collection at laminar, transitional, and 
fully-turbulent (untripped) conditions. Laminar and transitional data were collected at Mach 10 conditions. 
Laminar through fully-turbulent data were obtained at Mach 8 conditions. Perfect gas heating data are 
valuable for code validation because the heating is a result of conduction only and does not include any 
diffusion effects (Equation 1). The AEDC test provided an excellent opportunity for CFD code comparison 
to laminar through turbulent data in a single facility. 

LAURA and DPLR solutions were obtained for low (laminar) and high (turbulent) Reynolds numbers 
at AEDC facility conditions. The codes were run with perfect gas nitrogen and a fixed wall temperature. 
As with flight predictions, the Baldwin-Lomax solutions were run with transition specified at the stagnation 
point. The LAURA and DPLR laminar heating compares very well to data at a low Reynolds number 
(Figure 6a). Both codes are within 10% of the time-averaged heat transfer at every measurement location. 
The error bars on the data (+/- 10%) are shown for reference purposes and do not necessarily represent 
the actual experimental uncertainty. LAURA and DPLR Baldwin-Lomax turbulent heating at a higher 
Reynolds number match the data almost as well as the laminar case, especially on the leeside (Figure 6b). 
The DPLR SST model predicts slightly higher leeside flank heating than both Baldwin-Lomax predictions, 
as well as the data. The code-to-code differences using the Baldwin-Lomax model are considered to be within 
normal CFD uncertainties. The small differences between the DPLR Baldwin-Lomax and SST models are 
also considered insignificant. It should not be surprising that the Baldwin-Lomax algebraic model predicts 
blunt body turbulent heating as well as the two-equation SST model, given that the flowfield is completely 
attached and does not require advanced turbulence modeling. Overall, the good comparisons between the 
CFD codes and AEDC data gave confidence in using the Baldwin-Lomax model for MSL turbulent heating. 




(a) Laminar ( Mach = 10, Re^p = 0.93 X 10 6 , a = 16) 


(b) Turbulent ( Mach = = 24.8 X 10 6 ,a = 16) 


Figure 6. CFD Compared to AEDC Tunnel 9 Data 


A heating augmentation was present in the stagnation data from both AEDC and the NASA Langley 
Research Center’s 20-Inch Mach 6 Air Tunnel. 20 For example, Figure 6b shows computed stagnation heating 
more than 10% below the data. The origin of the heating discrepancy is unknown, but it may be related to 
turbulence or unsteady acoustic feedback from the boundary layer edge. The discrepancy between computed 
and measured stagnation heating in Figure 6b may indicate the beginnings of that phenomenon. The 
heating bump exhibited in Figure 6b and other data has been included in the MSL design environments as 
an uncertainty near the stagnation point. 
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C. CFD Code Comparison at Flight Conditions 

No Mars flight data exist that can be used to validate turbulent heating predictions at flight conditions. 
Pathfinder thermocouple data 37 indicate laminar conditions or possibly early onset of transition at the 
shoulder, but no definite signs of turbulence. Given the lack of Mars flight data, CFD code comparisons 
were the next best option to verify the Baldwin-Lomax model implementations at MSL flight calculations. 
The entry trajectory (07-25) that was used for this exercise is analogous to 08-TPS-02, but with a 5.8 km/s 
entry velocity instead of 5.9. Figure 7 shows a comparison of baseline LAURA and DPLR Baldwin-Lomax 
environments at times of peak heat flux, peak shear stress, and peak pressure on the 07-25 trajectory. At 
the leeside shoulder, the heat flux (Figure 7a) and total heat load (Figure 7d) results are nearly identical 
between codes. The discrepancies at other locations are considered to be within normal CFD modeling 
uncertainties. The shear stress agreement is similar to heat flux (Figure 7b) and the pressures match very 
well (Figure 7c). Based on the good agreement between LAURA and DPLR, no uncertainties were added 
to account for code-to-code Baldwin-Lomax model variability in the flight design environments. 




X (m) X (m) 

(a) Heat Flux (b) Shear Stress 


1 07-25 Nominal Trajectory (Baseline) | 



X (m) 



X (m) 


(c) Pressure 


(d) Total Heat Load 


Figure 7. LAURA and DPLR Baldwin-Lomax Solutions on the 07-25 Trajectory (No Uncertainties) 
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D. Aerothermodynamic Partials 

There are dozens of models in LAURA and DPLR needed to simulate the chemically and thermally non- 
equilibrium flowfield that is predicted to surround the MSL heatshield at flight conditions. These models are 
governed by numerous parameters and often have questionable applicability for Mars hypersonic conditions. 
Some phenomena that affect aerothermodynamics simply may be beyond state-of-the-art CFD capabilities. 
The expectation of turbulent conditions during flight adds yet another source of modeling uncertainty that 
pushes the limits of existing CFD models. In light of the numerous models that affect CFD predictions, an 
attempt was made to quantify some of their uncertainties and effects on the MSL aerothermodynamic design 
conditions. Previous attempts 38 have been made to quantify CFD uncertainty effects on Pathfinder heating 
calculations. The MSL project had neither the time nor the resources to perform a rigorous uncertainty 
analysis. Rather, a limited set of models and phenomena that were deemed important were considered. For 
the most part, conservative assumptions (e. g. a fully-turbulent heat pulse) were chosen up front to avoid 
having to estimate more detailed CFD uncertainty effects on the MSL heatshield design conditions. 

The MSL design heatshield environments were specified to be the baseline LAURA results using the 
stated model assumptions plus partials ( biases and uncertainties ) obtained through code validation and 
sensitivity analyses: 


( q,r,p) Des ign = ( q,T,p) L AURA + A(q,T,p) Bias + A(q, T,p) Uncert (4) 

Biases are corrections to the heat flux, shear stress, and pressure due to CFD model deficiency or inaccuracy; 
they are applied sequentially to the baseline environments: 


A(q,T,p) B ias 


d(q, T,p) w 
^ dXi 


(5) 


where represents those parameters that are expected to alter the environments. Uncertainties are those 
phenomena that may change the environments; they are treated as random independent occurrences and are 
combined using a Root Sum Square (RSS) approach: 


A (q,r,p) Uncert 


V ^ V dXi 


(6) 


Table 3 summarizes the partials and applicable parameter ranges. The Surface Catalysis bias reflects 
the use of the super- catalytic boundary condition, which forces the recombination of CO 2 and N 2 to their 
freestream mass fractions. In the absence of reliable data, the super-catalytic assumption is frequently used 26 
because it gives the highest surface chemistry contribution to heat flux (Equation 1), and thus is conserva- 
tive for TPS design. Figure 8 illustrates the conduction and diffusion contributions to super-catalytic heat 
flux on the 08-TPS-02 trajectory. At the time of peak heating, diffusion is predicted to contribute about 
two-thirds of the total heat flux at the leeside shoulder. The individual contributions are about equal near 
the stagnation point. Pathfinder thermocouple data suggests that, for those flight environments, the heating 
was 85-100% of super-catalytic CFD computations. 37 The Turbulent Augmentation bias (5%) reflects the 
small differences between CFD and AEDC turbulent data on the leeside flank (Figure 6b). The Distributed 
Roughness bias is discussed below. 

LAURA solutions were used to determine most of the uncertainties by individually varying the indepen- 
dent parameter and calculating the effect on the environments. Parameter ranges were taken from either 
separate analyses (e. g. Trim Angle- of- Attack and Shape Change ) or expert judgment (e. g. Turbulent 
Schmidt Number). Differences between the nominal and dispersed environments (heat flux, shear stress, 
and pressure) were used as the uncertainties at a given heatshield location. All CFD-based uncertainties 
were calculated at the time of peak heat flux and summed using the RSS method (Equation 6) before being 
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Table 3. Aerothermodynamic Partials 


Type 

Independent Parameter 

Baseline 

Minimum 

Maximum 


Surface Catalysis 

Super- Catalytic 

N/A 

N/A 

Bias 

Distributed Roughness 

Eqns. 7, 9 (k s = 0.6 mm) 

N/A 

N/A 


Turbulent Augmentation 

+5% 

N/A 

N/A 


Trajectory 

Nominal 

-3 a 

T3 a 


Turbulent Schmidt Number 

0.7 

0.5 

1.0 


Trim Angle-of- Attack 

Nominal a tr im 

(%trim “ 5 

(%trim T 5 


Surface Emissivity 

0.89 

0.5 

1.0 

Uncertainty 

Transition Location 

Xtran 0 

X-tran 0 

A/ran = -D/4 

Distributed Roughness 

0 

-10% 

+10% 


Stagnation Augmentation 

None 

None 

+50% 


Shape Change 

None 

None 

Deformed 


Parametric (q w ,r w ,p w ) 

0 

-(10,10,4)% 

+(10,10,4)% 


Other ( q w ,r w ,p w ) 

0 

-(5,5,2)% 

+(5,5,2)% 


applied uniformly along the trajectory. Shape Change refers to the effects of TPS recession and aeroshell 
deformation on the aerodynamic surface. CFD analysis was performed to determine the environments on a 
deformed shape with a maximum PICA recession of 0.8 in on the leeside shoulder. Turbulent heat fluxes on 
the deformed and recessed shape were within 0.5 % of each other, an amount that was considered negligible. 
The Stagnation Augmentation uncertainty refers to MSL experimental data 16 in which the stagnation point 
exhibited a heat flux bump at high Reynolds numbers for unknown reasons. This partial was only applied to 
those points near the stagnation point. The Parametric uncertainty was included to represent computational 
modeling inputs e. g. reaction rates, transport properties. Other represents an estimate of uncertainties 
that are not explicitly addressed. 

All CFD cases were run assuming a smooth 
aerodynamic surface representing the heatshield 
outer mold line. PICA will char and ablate 
when exposed to the expected heating in flight, 
resulting in a roughened surface. If the sur- 
face roughness elements extend far enough into 
the turbulent boundary layer, heat flux and 
shear stress will be augmented above smooth 
surface levels. The augmentation of turbu- 
lent heat flux and shear stress has been docu- 
mented many times in experiments . 39,40 The Dis- 
tributed Roughness bias in Table 3 captures the 
expected turbulent heat flux and shear stress 
augmentation caused by an ablated PICA sur- 
face. Roughness effects for MSL were estimated 
with curve fits of experimental data 39 that corre- 
late the turbulent shear stress bump factor due 
to roughness and the roughness Reynolds number 

{k+y. 


08-TPS-02 Nominal T rajectory at Peak q w , LAU R A B-L Solutions 



Figure 8. Conduction and Super-Catalytic Diffusion 
Contributions to 08-TPS-02 Peak Heat Flux 


^ = 1+ 0.9 (log 10 k + - 1), 10 < k + < 70 (7) 

t w ,o 


where 
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PwU r k, 

pw 


(8) 


k + = 


and k s is the equivalent sand grain A0 roughness height. The heat flux bump factor due to roughness is 
derived 41 from Equation 7: 


1+0.6 

Qw, o \r w ,o 


(9) 


Figure 9a shows the skin friction coefficient (i. e. shear stress) bump factor data used for MSL. The range 
where 10 < k + < 70 is termed the transitionally rough regime. Equation 7 is simply a fit of the middle line 
in Figure 9a. The 08-TPS-01a roughness Reynolds numbers as computed with the LAURA Baldwin-Lomax 
model are shown in Figure 9b. Peak leeside shoulder values for k + are in the transitionally rough regime 
and will lead to augmented shear stress and heat flux. The correlations given by Equations 7 and 9 were 
applied to the LAURA Baldwin-Lomax smooth- wall results as a bias (Table 3). Roughness was also applied 
as an uncertainty (Table 3) to cover data scatter in Figure 9a. 


1 08-TPS-01 a Nominal T rajectory, LAU R A B-L Solutions | 




(a) Correlation of Shear Stress Bump Factor vs. Roughness (b) LAURA Roughness Reynolds Number Distributions on the 
Reynolds Number 39 08-TPS-01a Trajectory ( k s = 0.6 mm) 


Figure 9. Distributed Roughness Effects 


The assumption for MSL is that the PICA equivalent sand grain roughness is equal to the physical 
roughness height ( k s = k in Equation 8). No post-flight measurements of PICA roughness exist where the 
material has been exposed to turbulent conditions. Stardust post-flight measurements (laminar) showed 
roughness heights on the order of 0.1 mm. 6 MSL post-test measurements of PICA samples exposed to 
laminar and turbulent conditions showed a maximum roughness height of 0.6 mm ; this value was used for 
the MSL heatshield environments. Figure 10 shows the bump factors along the 08-TPS-01a trajectory using 
the 0.6 mm roughness height. The leeside shoulder bump factors are highest (large & + ), with maximum 
shear stress and heat flux factors of 1.40 and 1.24, respectively. The bump factors were applied directly 
to the LAURA smooth wall results prior to adding uncertainties. Given the impact of PICA roughness on 
the heatshield design environments, MSL convened an independent advisory group to review the analytical 
methods presented here and to recommend additional testing. 6 MSL responded to those recommendations 
that could be accommodated prior to the original 2009 launch date. 
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Another mechanism for heat flux augmentation is protruding gap filler material between PICA tiles. TPS 
testing at some conditions showed that the gap filler material recessed more slowly than the surrounding 
PICA, causing ridges between tiles. 7 This mechanism was not specifically included in the design environ- 
ments since the likelihood of protrusions occurring is low 7 in high heating areas and the effects on heating 
are localized. 25 



X (m) 


(a) Shear Stress 


|08-TPS-01a Nominal Trajectory, LAURA B-L Solutions | 



Figure 10. Bump Factor Distributions due to Distributed Roughness on the 08-TPS-01a Trajectory 


Figure 11a shows the 08-TPS-02 baseline (smooth- wall, super-catalytic) peak environments with the 
individual biases and uncertainties. Maximum baseline heat flux at the leeside shoulder is 128 W/cm 2 
prior to application of the biases and uncertainties. The largest leeside shoulder heat flux partial is the 
Distributed Roughness bias. The same is true for shear stress (not shown). Trim Angle- of- Attack, Trajectory , 
and Turbulent Schmidt Number are among the largest uncertainty contributors. The total heat flux bias 
(Equation 5) and uncertainty (Equation 6) are shown in Figure lib. At the leeside shoulder, the total bias 
and uncertainty magnitudes are 32 W / cm 2 and 36 W/cm 2 , or 25% and 28% of the baseline heat flux. At the 
same location, the cumulative bias plus uncertainty (68 W/cm 2 ) is about 53% of the baseline. Note that the 
heat flux uncertainty associated with the stagnation region augmentation is not represented in Figure 11. 



X (m) 


1 08-TPS-02 Heatshield Peak Heat Flux Partials | 



(a) Individual 


(b) Total 


Figure 11. 08-TPS-02 Peak Heat Flux Biases and Uncertainties 
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E. Design Environments 


Figure 12 and Table 4 summarize the peak design environments and conditions at the peak locations, regard- 
less of when they are predicted to occur. Peak heat flux, shear stress, and pressure values correspond to the 
steep 08-TPS-02 trajectory. Design total heat load corresponds to the shallow 08-TPS-01a trajectory. The 
design values in Table 4 are simply the peak values regardless of trajectory. The Stagnation Augmentation 
uncertainty for heat flux is not shown in Figure 12a. Turbulent heat flux and shear stress are highest at 
nearly the same location on the leeside shoulder. At that location, the design values are 197 W/cm? for 
heat flux, 471 Pa for shear stress, 0.262 atm for pressure, and 5477 J /cm? for total heat load. Maximum 
design stagnation point pressure is 0.371 atm. The Leeside Flank location requires the most PICA thick- 
ness (0.94 in) because the underlying structural density (i. e. thermal mass) is reduced radially outward 
starting at that location. 6 The maximum heat flux location has more severe environments, but they are 
highly localized and are offset by lateral conduction relief in the shoulder. MSL’s design heating conditions 
far exceed what Pathfinder was estimated to have experienced 37 (118 W/cm? and 3800 J /cm? assuming 
laminar super-catalytic conditions). The shear stress and pressure conditions are also unprecedented for a 
Mars entry capsule. 



o 

X (m) 


(a) Symmetry Plane Heat Flux, Shear Stress, and Pressure 


08-TPS-02 Heatshield Design Peak Heat Flux Environments 



(b) Heat Flux 


1 08-TPS-02 Heatshield Design Peak Shear Stress Environments | 




08-TPS-02 Heatshield Design Peak Pressure Environments 


P w (atm) 


u 


0.371 

0.334 

0.297 

0.260 

0.223 

0.185 

0.148 

0.111 

0.074 

0.037 

0.000 



(c) Shear Stress 


(d) Pressure 


Figure 12. Heatshield Design Peak Environments 
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Table 4. Heatshield Design Environments 


Location 

Peak q w ( W j cm 2 ) 

Peak r w (Pa) 

Peak p w (atm) 

Qw ( J/cm 2 ) 

Maximum Heat Flux 

197 

444 

0.262 

5477 

Maximum Shear Stress 

178 

471 

0.227 

5054 

Maximum Pressure 

45 

10 

0.371 

2219 

Leeside Flank (X = -1.8 m) 

191 

372 

0.268 

5224 


Figure 13 shows the time-dependent design conditions at the maximum heat flux location on both design 
trajectories. The heat flux, shear stress, and pressure pulses are offset from one another, so peak values will 
not occur simultaneously. At all heatshield locations, the heat flux peak occurred first, followed by shear 
stress and pressure. The leeside shoulder is predicted to experience maximum heat flux and shear stress at 
essentially the same location and within a few seconds of one another, which will be particularly stressful for 
the PICA material. Simultaneously high heat flux and shear stress could lead to ablation and mechanical 
erosion of PICA. Previous entry capsules whose peak heat flux occurred under laminar conditions (e. g. 
Pathfinder) did not experience significant shear stress at the same location. 

Time histories like those shown in Figure 13 were provided for TPS modeling to predict the amount of 
mass blowing, surface recession, required thickness, etc. during flight. The TPS model includes material 
response biases and uncertainties analogous to those used for the aerothermodynamic environments. At the 
limiting location, 0.94 in of PICA are required to maintain a bondline temperature of 250° C, 6 including 
all thermal analysis biases and uncertainties. The PICA response for MSL is dominated by recession (0.73 
in maximum) due to the severe conditions and large model uncertainties. 6 The additional 0.31 in of PICA 
on the flight heatshield is considered insurance against potential risks in the predicted aerothermodynamics 
and TPS material response. 


08-TPS-01 a Design Environments: HS Leeside Shoulder (Max. q w ) 08-TPS-02 Design Environments: HS Leeside Shoulder (Max. q w ) 
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(a) 08-TPS-01a Trajectory 


(b) 08-TPS-02 Trajectory 


Figure 13. Heatshield Design Environments at the Maximum Heat Flux Location 
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Aerothermodynamic environments were also de- 
rived for the development of the MSL Entry, 
Descent, and Landing Instrumentation 9 (MEDLI) 
sensors. MEDLI consists of two suites of in- 
strumentation: MEDLI Integrated Sensor Plugs 

(MISP) temperature/recession sensors and MEDLI 
Entry Air Data System (MEADS) pressure sen- 
sors. The sensor distributions are shown in Fig- 
ure 14. The sensors are installed into PICA plugs 
that are subsequently flush-mounted to the flight 
heatshield. The data collected by MEDLI will be 
used to reconstruct surface heat flux (MISP) and 
capsule surface pressure and attitude (MEADS). 
The MISP locations cover a broad range of ex- 
pected heat flux environments, including near the 
maximum heat flux and shear stress on the lee- 
side shoulder (T2/T3). The MEADS locations 
are concentrated in the higher pressure and lower 
heat flux/shear stress region near the stagnation 
point and nose. Biases and uncertainties were 
developed for each MEDLI sensor location us- 
ing the method described previously. The de- 
sign conditions are summarized in Tables 5 and 
6 . 



Figure 14. MEDLI Sensor Locations 


Table 5. MISP Design Environments 


Location 

Peak q w (W / cm 2 ) 

Peak r w (Pa) 

Peak p w (atm) 

Qw ( J/cm 2 ) 

T1 

56 

8 

0.371 

2345 

T2/T3 

189 

376 

0.267 

5168 

T4 

52 

7 

0.370 

2092 

T5 

107 

124 

0.265 

3637 

T6 

184 

285 

0.264 

5028 

T7 

118 

101 

0.248 

3711 


Table 6. MEADS Design Environments 


Location 

Peak q w (W/ cm 2 ) 

Peak r w (Pa) 

Peak p w (atm) 

Q w (J/cm 2 ) 

pi 

52 

8 

0.369 

2117 

P2 

53 

6 

0.370 

2169 

P3 

72 

33 

0.362 

2812 

P4 

98 

83 

0.322 

3484 

P5 

107 

131 

0.245 

3591 

P6/P7 

95 

70 

0.295 

3263 
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TPS test facilities cannot simulate the time-dependent conditions that PICA will experience during flight. 
Rather, the facilities are limited to a small number of simultaneous surface conditions (heat flux, shear stress, 
and pressure), some of which may not be representative of flight. Given this limitation, the MSL TPS test 
program focused on testing at conditions at which one, or perhaps two, of the surface conditions were sim- 
ulated at the same time. 7 Priority was given to testing in shear conditions to simulate the simultaneously 
high heat flux and shear stress on the leeside shoulder. Figure 15 shows envelopes of peak conditions that 
were used to select facilities and test environments. The MSL program used these envelopes to determine 
test conditions for development and qualification of PICA (Figure 15a) and the MEDLI sensors (Figure 15b 
and Figure 15c). MISP development required testing near the peak values for both heat flux and shear stress 
at the same time. MEADS testing focused its testing in high pressure conditions with moderate heat flux 
and low shear stress. 


08-TPS-02 Heatshield Design Environments Envelope | 


08-TPS-02 MISP Design Environments Envelope 



(a) Peak Locations 



(b) MISP Locations 


08-TPS-02 MEADS Design Environments Envelope 



P w (atm) 


(c) MEADS Locations 


Figure 15. Heatshield Design Environments Envelopes 
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F. Aerothermodynamic Indicators 

The computation of aerothermodynamic environments via Navier-Stokes flowfield solutions consumes signif- 
icant time and resources. Any changes to the entry capsule mass properties, atmospheric density profile, 
launch dates, etc. feed into the design trajectory and require a repeat of the aerothermodynamic design 
process (Figure 4). In order to avoid having to run CFD for every candidate trajectory, a method was 
employed in which environments were estimated using curve fits of CFD results from a reference trajectory. 
This method follows previous work 42 in which CFD-based surface environments are fit to an equation that 
depends only on freestream density and velocity: 


q w ,r w ,p w = CpZV£ 


( 10 ) 


where the units are kg/m 3 for and m/s for V^. The coefficient (C) and exponents (m, n ) were derived 43 
to fit Equation 10 to time- varying design conditions, such as those shown in Figure 13. 


Figure 16 shows 08-TPS-01a design condi- 
tions and derived curve fits ( aerothermodynamic 
indicators ) at the maximum heat flux loca- 
tion. The curve fitting parameters are shown 
in Table 7. The curve fits were anchored to 
the peak CFD values by adjusting the C co- 
efficients. The benefit of these indicators is 
that they can be tracked in trajectory simula- 
tions to estimate aerothermodynamics on alter- 
native trajectories prior to (or instead of) run- 
ning CFD solutions. For instance, the effects 
of landing location or arrival date could be es- 
timated in a fraction of the time it would take 
to run the CFD solutions by applying the in- 
dicators to the new trajectory. Of course, 
the CFD solutions should be updated for TPS 
analysis once a final design trajectory is se- 
lected. 
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Figure 16. 08-TPS-01a Design Environments and An- 
chored Indicators at the Maximum Heat Flux Location 


Table 7. 08-TPS-01a Indicators at the Maximum Heat Flux Location 


q w ,T w ,p w = Cp™V™ 

C 

m 

n 

q w ( W/cm 2 ) 

0.224 x 10" 10 

0.879 

4.22 

t w (Pa) 

0.198 x 10 -8 

1.082 

3.95 

p w (atm) 

0.135 x 10 -5 

1.055 

2.26 
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G. Design Environments for a 2011 Launch 


Aerothermodynamic indicators were used to pre- 
view the impact of delaying MSL launch until 2011. 
Until late 2008, launch was scheduled for late 2009. 
For various technical and programmatic reasons, the 
decision was made to delay launch until the next 
opportunity in late 2011. The major effect of the 
delayed launch is an increased entry velocity from 
5.63 to 5.90 km/s. The actual magnitude of this 
velocity increase will vary with arrival date. Be- 
fore CFD analysis began for the 2011 launch oppor- 
tunity, aerothermodynamic indicators derived from 
the 08-TPS-01a/02 trajectories were used to esti- 
mate the new environments on an initial 2011 design 
trajectory. Figure 17 shows indicators derived from 
08-TPS-01a at the maximum heat flux location ap- 
plied to the 09-TPS-01 (shallow entry from a 2011 
launch) trajectory. The 09-TPS-01 CFD-based de- 
sign environments have since been completed using 
the process described previously. Figure 17 shows 
that the indicators accurately predicted the CFD- 
based environments. The peak indicators values are 
all within 10% of the CFD peaks. 
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Figure IT. 09-TPS-01 (2011 Launch) Indicator-Based 
and CFD-Based Design Environments at the Maximum 
Heat Flux Location 


Table 8 compares the design environments at critical locations for the 2009 and 2011 launches. The 
change in launch date increases turbulent heat flux, shear stress, and total heat load on the leeside shoulder. 
Conversely, pressure decreases at all locations. An initial TPS assessment 6 has shown that an additional 
0.07 in (1.02 in total) is needed to accommodate the increased heating. As new trajectory analyses are 
completed for the new launch date, the aerothermal team will revisit the environments analysis to assess the 
risks introduced and their impacts on the current PICA thickness. 


Table 8. Heatshield Design Environments for the 2009 / 2011 Launch Opportunities 


Location 

Peak q w ( W j cm 2 ) 

Peak r w (Pa) 

Peak p w (atm) 

Qw ( J/cm 2 ) 

Maximum Heat Flux 

197 / 226 

444 / 466 

0.262 / 0.243 

5477 / 6402 

Maximum Shear Stress 

178 / 203 

471 / 490 

0.227 / 0.209 

5054 / 5895 

Maximum Pressure 

45 / 48 

10/5 

0.371 / 0.332 

2219 / 2042 

Leeside Flank (X = -1.8 m) 

191 / 220 

372 / 392 

0.268 / 0.246 

5224 / 5976 


III. Summary and Conclusions 

Aerothermodynamic design conditions (surface heat flux, shear stress, pressure, and total heat load) 
were derived for the Mars Science Laboratory capsule heatshield for entry following the original 2009 launch. 
Navier-Stokes flowfield solutions obtained on steep and shallow design trajectories formed the basis of the 
environments. Chemical and thermal non-equilibrium models were used to simulate the hypersonic and high- 
temperature conditions that will envelope the entry capsule during Mars atmospheric entry. The expectation 
of boundary layer transition in flight, a first for Mars entry capsules, necessitated the prediction of fully- 
turbulent conditions for the design heat pulses. The flowfield solvers using the Baldwin-Lomax algebraic 
turbulence model agree well with fully-turbulent experimental heating data. Code-to-code discrepancies for 
flight environments predictions were within normal CFD uncertainties. Biases and uncertainties were added 
to the baseline CFD results to represent computational uncertainties or phenomena not included in the 
CFD models. A correction bias was included in the design conditions to account for heatshield TPS surface 
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roughness effects on turbulent shear stress and heat flux. At the maximum heat flux location, the heat flux 
total bias plus uncertainty is more than 50% of the baseline heating. 


The design environments were taken to be the maximum values from either the steep (maximum heat 
flux, shear stress, and pressure) or shallow (maximum total heat load) design trajectory. The peak design 
conditions regardless of location are: 197 W/cm 2 heat flux, 471 Pa shear stress, 0.371 atm pressure, and 
5477 J /cm 2 total heat load. Location-specific conditions were derived as a function of time and used in 
modeling and testing of the heatshield TPS material. A method of estimating aerothermodynamic conditions 
on alternative trajectories was developed using curve fits of existing CFD-based environments. Finally, 
delaying MSL launch until 2011 will increase design peak heat flux and total heat load mostly due to a 
higher atmospheric entry velocity. 
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